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ABSTRACT 

The high-frequency-peaked BL Lacertae object Markarian 180 (Mrk 180) was observed to have an 
optical outburst in 2006 March, triggering a Target of Opportunity observation with the MAGIC 
telescope. The source was observed for 12.4 hr, and very high energy 7-ray emission was detected 
with a significance of 5.5 a. An integral flux above 200 GeV of (2.3 ± 0.7) x 10~^^cm~^s~^ was 
measured, corresponding to 11% of the Crab Nebula flux. A rather soft spectrum with a photon 
index of —3.3 ± 0.7 has been determined. No significant flux variation was found. 
Subject headings: BL Lacertae objects: individual (Markarian 180) — gamma rays: observations 



1. INTRODUCTION 

The search for very high energy (VHE, defined as 
E > 100 GeV) 7-ray emission from active galactic 
nuclei (AGNs) is one of the major goals for ground- 
based 7-ray astronomy. New detections open the pos- 
sibility of a phenomenological study of the physics in- 
side the relativistic jets in AGNs, in particular, to un- 
derstand both the origin of the VHE 7-rays as well 
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as the correlations between photons of different energy 
ranges (from radio to VHE). The number of reported 
VHE 7-ray-emitting AGNs has been slowly increasing 
and is currently 12 (2006 June). Six of them hav e 
been seen by MAGIC: Mrk 421 l|Albert e t al.l l2006d) 
Mrk 501 (Mazin 200i lES 1 9 59-^6 50 (Albert et all 
I2006a|). 1 ES2344+514 (IMazin I I200I . 1ES1218+304 
^^ertet al. 2006l3), and PG1553-hll3 (jXlbcrt et al] 
!2006d). 

The known VHE 7-ray-emitting AGNs are variable 
in flux in all wave bands. Correlations b etween X-ray 
and 7-ray emission have been found (e.g., iFossati et alJ 
(2004)), although the relationship has proven to be 
rather complicated, with 7-ray fla res also beirig de- 
tected in the absence of X-ray flares (iHolder et al.ll2003l: 
Krawczvnski et al1l2003() and vice versa (jRebillot et alJ 
2003ft . The optical- TeV correlation has yet to be stud- 
ied, but the op tical- GeV correlations seen in 3C 279 
()Hartman et al .' 2001) suggest that at least in some 
sources, such correlations exist. Using this as a guide- 
line, the MAGIC collaboration has been performing Tar- 
get of Opportunity observations whenever alerted that 
sources were in a high flux state in the optical and/or 
X-ray band. 

The AGN Mrk 180 (lES 1133-1-704) is a well-known 
high-frequency-p eaked BL Lac (H BL) object at a red- 
shift of z = 0.045 ()Falco et al.ll999j) . The spectral energy 
distribution (SED) of HBL objects exhibits a generic two- 
bump structure: one peak with a maximum in the X-ray 



2 



band and the other peak located in the GeV-TeV band. 
The radiation is produced in a highly beamed plasma 
jet, which is almost aligned with the observer's line of 
sight. A double-peaked SED is normally attributed to 
a population of rclativistic electrons, where one peak 
is due to synchrotron emission in the magnetic field of 
the jet and where the second peak is caused by inverse 
Compton (IC) scattering of low-energy photons. The 
low-energy photons can be external to the jet (e xternal 
Compton scattering : lDermerfc Schlickeise r' (1993)) or are 
produced within the jet via synchrotron radiation (syn- 
chrotr on self-Compton [SSC] scattering; iMarasch i et al.l 
l)1992l) \ Models Isased on the acceleration of hadrons 
can als o sufficiently describe the observed SEDs and light 
curves l|Mannheim et alJ[T99lt iMassaro et alJl2004 1. AU 
of the up to now known AGNs with strong GeV/TeV 
7-emission belong to the HBL object class. 

In many cases, the second peak of the SED is not 
observable because of low sensitivity above a few times 
100 MeV of satellite-borne detectors and a too high en- 
ergy threshold of ground-based 7-ray detect ors. In case 
of Mrk 180, HE GRA ( Aharonia n^ral]l200|) and Whip- 
ple ijHoran et alj.2004) observed this object but were only 
able to derive flux upper limits, and EGRET did not de - 
tect the source ijFichtel et al.ll994 >'Hart man et alll999|) . 

In the optical, Mrk 180 is characterized by a bright 
host galaxy R — 14.17 ± 0.02 mag and a much fainter 
(varia ble) core R = 15.79 ± 0.02 mag ijNilsson et al.l 
1200.1(1 . The variations in total optical brightness are 
therefore small. The source is observed regularly as part 
of the Tuorla Observatory Blazar monitoring program^'^ 
with the Tuorla 1 m telescope and the KVA 35 cm tele- 
scope^'*. To determine the core flux, we had to subtract 
the flux of the host galaxy and a nearby star within the 
5" aperture radius (together R=14.96 mag, or 3.2 mJy; 
Nilsson, private communication). 

After Mrk 180 had been detec ted in the X-rays by 
HEAO-1 l|Mufson fc Hutted 1198 l(l . it was observed by 
various satellites with fluxes ranging from 6.3 to 22 x 
10~*^ ergs cm~^ s~^ ('Do nato et al.l l2?fflltl at around 
1 keV and from 5.0 to 9.8 x 10~^^ ergs cm~^ s~ * in 
the 2-10 keV band l|Perlman et all 120051 iDonato et all 
linnS). The source is monitored by the all-sky monitor 
(ASM) on board the Rossi X-Ray Timing Explorer. 

In this Letter we present the first detection of VHE 
7-ray emission from Mrk 180. 

2. OBSERVATIONS AND DATA ANALYSIS 

The MAGIC telescope (M ajor Atmosphe ric Gamma 
Imaging Cherenkov telescope; Toren j ()2004D ^ is located 
on the Canary Island La Palma (2200 m above sea level, 
28?45' north, 17?54' west). MAGIC is currently the 
largest imaging atmospheric Cherenkov telescop e with 
a 17 m -diameter tessellated reflector dish ( Cortin a~et al.l 
|2P05). The 3? 5 field of view camera comprises 576 pho- 
tomultipliers (PMTs) with enhanced quantum efficiency. 
The accessible energy range spans from 50-60 GeV (trig- 
ger threshold at small zenith angles) up to tens of TeV. 
The 7 point-spread function is about 0?1. 

The observation of Mrk 180 was triggered by a bright- 
ening of the source in the optical on 2006 March 23. The 




Fig. 1. — fl-band light curve of Mrk 180 extending from 
2003 January to the end of 2006 March as measured by the 
KVA telescope; 5 mJy is equivalent to 14.47 mag. 



E 

=" 2,5 



c 

— 



3; 5.5 
E 

■55.3 





1 


MAGIC(>200GeV) 

1 1 f 1: 


-T \ ASM - 

11,11 , , * . . 


i 

i 


i 


optical 



53824 

Time [MJD] 



23 See 

24 See 



htt2i//usersMitu^fi/kan 



http://tur3.tur. iac.es/l 



Fig. 2.— Light curve of Mrk 180 for MJD = 53815- 
53825 (March 21-31). Upper panel: VHE 7-rays measured 
by MAGIC above 200 GeV. Middle panel: ASM count rate. 
Lower panel: Optical flux measured by KVA. 



alert was given as the core fiux increased by 50% from its 
quiescent level value as determined from over 3 years of 
data recording, shown in Fig.^ During the MAGIC ob- 
servation, optical follow-up observations were performed 
with KVA. The simultaneous MAGIC, ASM, and KVA 
light curve is shown in Fig.|21 Around this time, Mrk 180 
was also observed as part of the AGN-monitoring pro- 
gram by the University of Michigan Radio Observatory 
(UMRAO). No evidence of flaring was found between 
2006 January and April. 

Mrk 180 was observed in 2006 during 8 nights (from 
March 23 to 31), for a total of 12.4 hr, at zenith an- 
gles ranging from 39° to 44°. The observations were per- 
formed in the so-called wobble mode (Damn et al. 1997), 
in which the telescope is pointed alternatively for 20 min- 
utes to two opposite sky positions at 0?4 off the source. 
Runs with unusual trigger rates due to detector problems 
or adverse atmospheric conditions were rejected. The to- 
tal observation time was thus reduced to 11.1 hr. 

The data were analyzed using the standard analy- 
sis and calibration programs for the MAGIC tele scope 
llBretzll2005t IWagner et al]l2005HGaug et al.ll2005l) . For 
the 7/hadron shower separation a multidimensional clas- 
siflcation t echnique based on t he Random Forest method 
was used l|Breimanl l200lt iBock et al.ll2004|) . The cuts 
for the 7/hadron separation were trained using a frac- 
tion of randomly chosen data to represent the back- 
ground (hadrons) and were compared to Monte Carlo 
(MC) 7 events with the same zenith angle d i stribu - 
tion (CORSIKA ver. 6.023; Kna nn HeTkl (|2004); 
iMaiumdar et alJ l|2005l) ). The cuts were then chosen 
such that the overall cut efficiency for MC 7 events was 
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Fig. 3.— The ^^-distribution for the ON -source data [filled 
circles) and normalized OFF-source events {gray histogram) 
for Mrk 180. The vertical dotted line indicates the 9^ cut used 
to determine excess events. The total excess of 165 events 
corresponds to a significance of 5.5 cr. 



> 1" r 
1 10-'° g- 

I 10' 

§ 10" 
z 

10" 



10"' 



dN/dE = f„*(E/0.3TeV) 




^f„ = (0.45 ±0.18)* 10 
i a = (-3.25 ±0.66) 



cm-^ s ' TeV' 



100 



1000 



E [GeV] 



Fig. 4. — Differential energy spectrum of Mrk 180. Full 
circles: Spectrum measured by MAGIC. Open circles: De- 
absorbed energy spectrum (see text). The horizontal bars in- 
dicate the size of each energy bin. The black line represents 
a power-law fit to the measured spectrum. The fit parame- 
ters are listed in the figure. For comparison, th e Crab Nebula 
energ y spectrum as derived from MAGIC data IWagner et alj 
120051) is shown {dotted line). 



about 50%. The number of excess events is determined 
as the difference between the source and background re- 
gion in the 9^ distributions, with 9 being the angular 
distance between the source position in the sky and the 
reconstructed arrival position of the air shower. The 
latter position is determined for each shower image by 
means of t he so-called DISP method (Fomin et al. 1994; 
iLessard'et al. 2001; Domingo-Santamarfa et al. 2005), 
using image shape parameters (.Hillasiiil985il . In order 
to determine the background distribution, three back- 
ground regions of the same size are chosen symmetrically 
to the on-source region with respect to the camera cen- 
ter. A final cut of 6*^ < 0.024 is applied to determine 
the significance of the signal and the number of excess 
events. The energy of the 7-ray candidates was also esti- 
mated using the Random Forest technique. The applied 
cuts were chosen to be looser than the ones in Fig. El in 
order to gain statistics on the 7-ray candidates. Due to 
the large zenith angle of the observation, the correspond- 
ing energy threshold (defined as the peak in the differ- 
ential energy distribution of the MC 7 events) after cuts 
was about 200 GeV. Effects on the spectrum determina- 
tion, introduced by the limited energy resolution of the 
detector, we re corrected using the "unfolding" methods 
according to lAnvkevev et aTl l)1991[l . 

3. RESULTS 

The distribution of 0^-values after cuts is shown in 
Fig. 131 The signal of 165 events over 605.2 normalized 
background events corresponds to a 5.5 a excess using 
equation (17) in iLi fc Mai l|198^ . The shape of the ex- 
cess is consistent with a pointlike source seen by MAGIC. 

No evidence of flux variability was found. The fit to the 
nightly integrated flux is consistent with a constant emis- 
sion (x^ = 7.1, 6 degrees of freedom). Fig. |21 shows the 
VHE light curve together with the ASM daily averages 
and the i?-band flux data. The X-ray flux of the source is 
generally below the ASM sensitivity, but on March 25 a 
3 a excess was observed, which suggests that the source 
was also active in X-rays. The optical flux reached its 
maximum on the night MAGIC started the observations 
(March 23) and began to decrease afterward. 

The measured energy spectrum of Mrk 180 is shown 



in Fig. 21 Assuming a power-law spectrum, we obtained 
the following parameterization: 

- = (4.5±1.8)xl0-"(^j TW^ 

The observed integral flux above 200 GeV is F{E > 
200 GeV) = (2.25 ±0.69) x lO^^cm-^s-i, which corre- 
sponds to 1.27 X 10~^^ ergscm^^s^\ or 11% of the Crab 
Nebula flux measured by MAGIC. The errors are statisti- 
cal only. We estimate the systematic errors to be around 
50% for the absolute flux level and 0.2 for the spectral 
index. The large systematic flux error is a consequence 
of the steep slope. The main error contributions are the 
uncertainty in the atmospheric transmission, the reflec- 
tivity and focusing uncertainty of the mirrors and light 
catchers in front of the PMTs, and the uncertainty in 
the effective quantum efficiency of the PMT and in the 
photoelectron-to-signal conversion. A second indepen- 
dent analysis gave results in very good agreement with 
the quoted numbers. 

The VHE 7-rays from Mrk 180 are partially ab- 
sorbed by the low-energy photons of the evolving 
extrag alactic background light (EBL; see INikishqvl 
llT96l: IGould Schroder (196fih: ISted^er et alJ l|199^ 
iHauser fc pTekl 1^20011) ). Given the redshift z=0.045 of 
Mrk 180, the effect is small for photons with energies 
below 1 TeV. We calculate the optical depth and the 
resulting attenuation of the VHE 7-rays from Mrk 180 
using the number dens ity of the evolving E BL provided 
by the best-fit model of Knei ske et al.l l)2004j) . This state- 
of-the-art model is consistent with the recently derived 
upper limits o n the EBL inf erred from arguments about 
AGN spectra llAharonian et'al]l2006D . The attenuation 
was det ermined by the numerical i ntegration of equation 
(2) from' Dwek fc Krennrichll|2005D . The de-absorbed en- 
ergy spectrum of Mrk 180 is also shown in Fig. 0] {open 
circles) . A fit with a simple power-law to the de-absorbed 
spectrum reveals a slope with a = —2.8 ± 0.7. 

4. DISCUSSION AND CONCLUSIONS 

In this Letter we have presented the first detection of 
VHE 7-ray emission from Mrk 180. The discovery was 
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Fig. 5.— SED of Mrk 180. Simultaneous data (UMRAO, 
KVA, ASM, and MAGIC) are represented as black circles. 
The gray circles rep resent historical data (Giomrni ct al. 2002; 
IPerlman et al.l200a. The arrows de note the upper limits from 
ASM. EGRET IFichtel et al.l 11991). Whipple and HEGRA. 
The solid line is fro ^Sostama^^^^^hisellinil I2002l') and the 
dashed line is from lFossat^^am2000l) (see textV 



triggered by an optical flare, but no significant varia- 
tions in the VHE regime were found. The short observa- 
tion period and the small signal prevent us from carrying 
out detailed studies. It is therefore impossible to judge 
whether the detected VHE flux level represents a flaring 
or a quiescent state of the AGN. 

Earlier observations by other experiments have only set 
upper limits on the VHE 7-ray emission from Mrk 180. 
HEGRA observed the source for 9.8 hr, resulting in an 
upper hmit (99% c.l.) at 1.5 TeV of 12% of the Crab 
Nebula flux f Aharonian et alJl2004|) . Whipple observed 
Mrk 180 during three observation periods for a total of 
26.8 hr, resulting in a flux up per limit at 300 GeV of 
10.8% of the Crab Nebula flux l|Horan et al.ll200^ . Both 
upper limits are above the flux presented in this Letter. 

Fig. El shows the SED of Mrk 180: the radio, opti- 
cal and X-ray data are from UMRAO, KVA, ASM, and 



the NED database, and the VHE data are from this 
analysis and various upper limits. Simultaneous data 
are noted in black, while historical data are noted in 
gr ay. The predicted flu x by the phenomenological model 
of iFossati et al.l ()200 n'l (dashed line) is too high, as al- 
ready emphasized by the non detection measurements 
l|Aharonia n et alJl2004t i Horan et afll2004D. A more de- 
tailed SSC model from iCostamante fc Ghisellinl l)2002f) 



[solid line) seems to describe the data more closely. 

The rather steep slope of the VHE spectrum suggests 
an IC peak position well below 200 GeV, while the non 
detection by EGRET gives a lower limit of ^ 1 GeV 
for the peak position. This is in agreement with the SSC 
model of Costamantc & Ghiscllini (2002) suggesting that 
the IC peak position is at around 10 GeV. The overall 
IC luminosity, however, is underestimated in this model: 
the observed integral flux above 300 GeV is a factor 30 
larger than predicted. An explanation for this discrep- 
ancy could be the model's underlying assumption of a 
quiescent-state synchrotron spectrum to obtain the IC 
flux. This could indeed suggest that our measurement 
was made during a high state. 

The discovery of VHE emission from Mrk 180 during 
an optical outburst makes it very tempting to specu- 
late about the connection between optical activity and 
increased VHE emission. Since Mrk 180 has not been 
observed prior to the outburst with MAGIC and since 
the upper limits from other experiments are above the 
observed flux level, further observations are needed. 
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